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We propose a method which can effectively stabilize fixed points in the classical and quantum 
dynamics of a phase-sensitive chaotic system with feedback. It is based on feeding back a selected 
quantum sub-ensemble whose phase and amplitude stabilize the otherwise chaotic dynamics. Al- 
though the method is rather general, we apply it to realizations of the inherently chaotic Ikeda map. 
One suggested realization involves the Mach-Zender interferometer with Kerr nonlinearity. Another 
realization involves a trapped ion interacting with làser fields. 
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Dctcrministic (Hamiltonian) chaos in a classical sys- 
tem manifests itself by the extreme (exponential) sen- 
sitivity of the evolution to initial conditions, thus mak- 
ing long-time predictions on the dynamical observables of 
the system practically impossible Q]. Yet, it has been re- 
cently recognized [|]|| that it is possible to control chaos: 
by applying specifically designed time-dependent pertur- 
bations, regular periòdic orbits can be stabilized within 
chaotic phase-space domains [|]|| ■ The quantum dynam- 
ics underlying classical chaos has been extensively investi- 
gated for some time now [0 . Of particular relevance here 
is the finding J5|] that quantum spread can cause a qua- 
siclassical state initially confined within regular phase- 
space regions to eventually access chaotic regions. It is 
therefore natural to ask: what is the quantum analog of 
classical chaos control and how can we deal with chaotic 
behavior induced by quantum spread? 

It is our purpose here to propose the quantum con- 
trol of deterministically chaotic two-dimensional systems, 
whose evolution is governed by their phase (mod 27r), 
and involves phase-dependent feedback. The proposed 
control is based on two elements. First, creation of phase 
correlations between the stabilized system and a similar 
system acting as a "stabilizer" . If the stabilized and sta- 
bilizcr systems are classically correlated and are not en- 
tangled, then their phase correlations allow us to enforco 
regular evolution (stabilization) of an inherently-chaotic 
stabilized system, by feeding it with the output of the 
stabilizer system. However, whenever these systems are 
appreciably entangled, this stabilization may not suffice. 
If we merely ignore (trace out) the stabilizer states, and 
feed back the resulting statistical mixture of stabilized- 
system states, then stabilization may eventually fail, be- 
cause one of the states in the mixture may venture into 
a chaotic region of phase space (see below). This calls 
for the second element of proposed control: a quantum 
measurement of the entangled stabilized and stabilizer 
systems, which seleets a particular state of the stabilized 



system and feeds it back into this system. The chaotic 
trend due to entanglement-induced spread is thus sup- 
pressed by suitable projections which localize the system 
in regular phase-space regions. Such manipulation of the 
evolution, i.e., selection of a certain state correspond- 
ing to the desired outeome of the measurement, is the 
essence of the conditional measurement (CM) approach 
to quantum state control || . The price one has to pay for 
guiding the evolution by CMs is that one must perform 
as many trials as implied by the success probability of the 
required CM until it is accomplished (see below). How- 
ever, we stress that this price is unavoidable, since there 
is no other remedy for the chaotic effects of entanglement. 
Since the quantum version of our scheme is non-unitary, 
it does not correspond to traditional control theory, in 
which a force driving the system can be identified. This 
essential difference stems from our objective, namely, to 
combat the unwanted effects of entanglement. We finally 
note that all of the above considerations apply on time- 
scales much shorter than dissipation (decoherence) time. 

Although the outlined proposal is quite general, we 
shall apply it here to systems governed by the Ikeda 
map (7,p. To this end, we shall analyze the Ikeda 
map quantum-mechanically, focusing on the quasiclas- 
sical limit, thereby extending its existing classical anal- 
ysis. Two diffcrcnt realizations will be considered: (i) a 
Mach-Zender interferometer with Kerr nonlinearity; (ii) 
a trapped ion interacting with traveling- and standing- 
wave fields. 

We first consider the following realization of the Ikeda 
map, which differs from its Standard realization in a ring 
cavity . The bàsic block of the proposed realization is a 
Mach-Zender interferometer (MZI) with a Kerr-nonlinear 
element in one arm, such that the output of one of the 
ports of this interferometer is fed back into the input 
port, along with a fixed input field (see block I in Fig. |^). 
The input-output operator transformation for this MZI 
is 
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FIG. 1. Schematic representation of the stabilized (Block I) 
and stabilizer (Block II) systems: two coupled Mach-Zender 
interferometers with Kerr nonlinearities (indicated by K) and 
feedback. The input ports are a and b. The 45° thin lines are 
50-50 beam splitters, while the 45° thick lines are perfectly 
reflecting mirrors. D denotes a homodyne detector. In the 
case of the single MZI the output / is fed back into port a 
(via the dashed-line beam splitter), while in the case of the 
total (system + stabilizer) setup the output r is fed back into 
a (via the thin-line beam splitter) . 
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where 0\ = (fi+Ki(b+iày (b+ia) is a phase-shift operator, 
with a c-number (refractive) shift ipi = X\h an d a Kerr 
shift which is proportional to the Kerr coefhcient n\ = 
Xj h/2 and to the field intensity d) d ||. 

The classical counterpart of this transformation, ob- 
tained by replacing all operators by c-numbers and set- 
ting 6 = (vacuum input in port b) can be used to de- 
scribe repeated feedback of the output into the input port 
a. If we choose to feed back the output / into port a to- 
gether with fixed external input (Fig.0), we obtain the 
following iterative map 



a in + ie i9l > /2 cos(%/2)a j 



(2) 



where j = 0, 1, 2, . . . is the iteration number, Oi n repre- 
sents the fixed input, and the 1/V2 factor is imposed by 
the 50-50 beam splitter at port a. 

This is a generalization of the Ikeda map in the limit 



of small Kerr shifts 



< <P, 



(3) 



with ài n = Oin/%/2, in that the feedback parameter R is 
now replaced by cos(0y/2)/·\/2, i-e., it varies with j. Its 



dynamics becomes chaotic whenever the c-number pa- 
rameter ài n /[l — | cos(0ij/2)|] is large enough 0, cor- 
responding to excessive feedback, and <p, k ^ [see 
Fig.l(a)]. 

We intend to show how the chaotic dynamics of this 
device can be made regular (stabilized) by connecting it 
to a similar block (block II in Fig. |ï|) via an interface. 
The second nonlinear MZI (block II) will then be shown 
to act as a "stabilizer" . The corresponding output-input 
operator transformations for the second block are 



= ie 
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where 6>n = ipu + Kn(f + iè)^(f + ie). 

In order to achieve stabilization we impose the condi- 
tions 
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(«i-kii) =Í(x[ 3) h-X y u>hi)=0. 
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(5a) 
(5b) 



At the level of classical wave òptics, when the opera- 
tors are replaced by c-numbers, these conditions ensure 
a fixed, j-independent (intensity-independent) feedback 
factor for the r-mode, namely, R — cos[(0t — 6*n)/2]/v / 2 = 
sm 5 /y/2. By choosing R < 0.1 we cause the dynamics 
to be stable [Fig. p|(b)] at the classical level. This re- 
sult implies that, classically, stabilization is achievable 
by additional rotation and phase shift of the two-mode 
basis, forcing the system from the unstable into the sta- 
ble domain. These additional rotation and phase shift are 
caused by an external "device" (the stabilizer) which, al- 
though similar to the controlled system, is not part of its 
dynamics, but rather the analog of the Ramsey apparatus 
which rotates the states of a two- level system M. If the 
controlled system is near the edge of the stability domain 
[cos(#i/2) just slightly larger than 0.14], then the required 
additional rotation and phase shift Ou/2 can be small and 
amount to a weak perturbation of the parameters, as in 
certain existing control methods ||. What singles out 
our classical stabilization method is that it exploits the 
interference of two correlated (coherent) evolutions. 

The more intricate part of our proposal is concerned 
with correcting for the unavoidable effects of quantum 
spread and entanglement, which cause the breakdown 
of the quantum-classical correspondence in phase space, 
when the Kerr effect is significant. Upon inverting Eq. (^) 
and using Eqs. (||), it can be readily verified that the 
output states corresponding to the modes r and q will in 
general be entangled, because of the cross-terms à^b+h.c. 
and thcir equivalent pè+h.c. in the phase-shift operators 

We are primarily interested in quasiclassical coherent- 
state inputs, which can be readily compared to their clas- 
sical counterparts 
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FIG. 2. (a) Chaotic dynamics of the íirst block. The 
"phase-space" picture [Re(a,j)-Im(aj)] is plotted for 10 iter- 
ations. The parameters used are tp = 7r/2 + tpi with ipi — 0.4, 
aín = 5.0 and k = 0.1. (b) Regular (stabilized) dynamics of 
the stabilized system enforced by the stabilizer system. Here 
the modulus of a-, is plotted versus j for the same vàlues of ai n 
and k as in (a), but with cp = n + (ipi + ipn) /2 = n/2 + ipi + S 
and R = 0.1. The fixed point (5.53,0.15) is reached after 
few iterations. Inset: Phase-plane distribution of a two-mode 
entangled superposition of coherent states for a fractional 
revival. Sòlid circles are "Schròdinger kittens" in one of 
the modes, with mean amplitude acosí 1 , dashed circles are 
their counterparts in the other mode, with mean amplitude 
asin«ï>. The dotted circle demarcates the boundary between 
the chaotic-domain and the regular-domain amplitudes. 



|in) = |a) a |0>6 = e 



_ P -M72 



n=0 



f\ra 



-|0>a|0)l 



(6) 



We then characterize our input density matrix by using 
the Q-function quasiprobability distribution, which is di- 
agonal in the coherent-state basis, that is, 
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Pf 



d 2 ad 2 0P(a,0)\a,0) a , b (a,0\ , 



(7) 



with P(a, 0) ~ g(a — a- ln , — 0), g being a narrow Gaus- 
sian. The resulting output, in general, will then be the 
entangled state 
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Pf 



d 2 ad 2 0P(a, /?)|out) r ,ç(out| 



(8) 



where the state |out) r . ç is in general an entangled state 
of field modes r and q. 

In order to achieve the same stabilization as in the 
classical limit, our goal is to project out at each iteration 
a coherent state, and feed it back from port r. It may be 
difficult to devise such a measurement always, since the 
entangled state (Jsj) may have a complicated phase-space 
distribution. We therefore choose to impose an additional 
condition on the Kerr shift, corresponding to a fractional 
revival of the Kerr shifted Fock states \n) [Tet] . The con- 
dition for a fractional revival is k = n\ + Kn = 7r(m//), 
where m and / are any integers. Then, for sufficiently 
large \a\ (\a\ ^> 1), state @ corresponds to an entangle- 
ment of a finite number of distinct coherent state compo- 
nents ( "Schròdinger kittens" ) in modes r and q with dif- 
ferent mean phases and amplitudes. A phase-plane plot 
of such a state reveals the origin of entanglement-induccd 
chaotic behavior [Fig. ||(b)-inset] : certain components of 
the fed-back mixture may be in the chaotic domain. Such 
a state allows us to project out the desired coherent-state 
component in r (localized in the regular domain), by a 
homodyne measurement of the g-component to which it 
is correlated. We illustrate this process by imposing the 
simplest (fractional) revival condition pï| 



K = Kj ■+■ Kll = 



1 n even 
i n odd 



(9) 



This condition is realizable for Kerr-shifts of photons in- 
teracting with àtoms in high-Q cavities ||, and, in a more 
straightforward fashion, for trapped ions interacting with 
làsers (see below) . Under this condition, we can show by 
a series of manipulations that Eq. (0) reduces to 
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|out) r>9 = -j= 



e -i7r/4| ae i* cos $) ç | - iae ^ S i n $) 



+e i7r / 4 | -io« i *am*> 9 |a!e i *cos$) r | , (10) 

where $ = {ip\ +(pn)/2. This entanglement of two-mode 
Schròdinger-cat states implies that we have 50% prob- 
ability of detecting the coherent state | — iae 1 * sin 
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Such detection is known to be possible by balanced homo- 



dyning, which mixes 
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local oscillator M. After a successful CM, we perforin a 
coherent mixing of the projected-out coherent state with 
the fixed input via a 50%-50% beam splitter (Fig. [[]) 
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finally feeding the resulting state [Eq. (|ll|) 1 back into the 
a-port. This result reduces [upon using the phase-space 
distribution (|§|)] to its classical counterpart Eq. (Q) for 
\a\ 3> 1. The above analysis demonstrates the possi- 
bility of achieving complete classical-like stabilization by 
an appropriate CM, thus overcoming the effects of strong 
quantum cntanglement. As stressed above, the price we 
have to pay is that we have to perform as many trials as 
implied by the success probability of the CM. 

The outlined scheme can be fully implemented in an 
ion trap, where the ion interacts with both traveling and 
standing-wave fields [jïï) . A two-mode classical traveling- 
wave pulse of area 9 can perform an adiabatic rotation of 
two orthogonal motional modes of the ion by the angle 
9: this is achievable via coupling to two orthogonally po- 
larized internal states, or by alternative schemes, based 
on measurements of the internal states |ïï|. The rota- 
tion of the basis of two motional modes thcn mimics the 
effect of beam splitters in the MZIs in Fig. [|. Kerr shift 
operators of the motional states nà^à are emulated by 
an off-resonant standing-wave field, whose k — Í7 2 r/A, 
Í2 being its Rabi frequency, r its duration, and A the 
detuning. Schròdinger-cat motional states have already 
been prepared in ion traps p2[ . 

To conclude, we have presented a scheme which en- 
forces regular dynamical evolution towards a stable fixed 
point in an inherently chaotic system by interference 
(phase correlations) with a "stabilizer" at the level of 
classical (wave) òptics. At the level of quantum descrip- 
tion, the present scheme is the first attempt to remedy 
the unwanted chaotic effects of quantum spread induced 
by the system evolution and its entanglement with a "sta- 
bilizer" device. This is achieved by a conditional (ho- 
modyne) measurement of the stabilizer output, which 
suppresses entanglement-induced chaos by projecting out 
a quasiclassical coherent state localized in the regular 
phase-space region. Conditions for such stabilization 
have been presented for the quasiclassical analog of the 
Ikeda map [Eqs. (|), @-(^)]. Although the present 
scheme is quite distinct from existing schemes, it is in 
the vein of those classical schemes of chaos control in 
which the phase space is drastically changed, either by 
an external (driving) force (the Pyragas method and its 
extensions) or by spectral filtering, for the purpose of 
stabilization ||. 
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